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xLiCl-(l  -x)(0.6Li20-0.4P205)  systems  with  x  =  0.1, 0.15, 0.2, 0.25,  have  been  prepared  using  melt  quench¬ 
ing  method  and  their  ionic  conductivity  was  characterized  by  impedance  spectroscopy.  Molecular 
dynamics  (MD)  simulations  for  the  same  systems  have  been  performed  with  an  optimized  potential,  fitted 
to  match  bond  lengths,  coordination  numbers  and  ionic  conductivity.  Based  on  the  equilibrated  config¬ 
urations  of  these  MD  simulations,  ion  transport  pathways  are  modelled  in  detail  by  the  bond  valence 
approach  to  clarify  the  influence  of  the  halide  dopant  concentration  on  the  glass  structure  and  its  conse¬ 
quence  for  Li  ion  mobility.  Features  of  the  consequential  ion  transport  pathway  models  (such  as  volume 
fraction  and  local  dimensionality  of  the  percolating  pathway)  are  compared  to  pathway  models  for  related 
glassy  solid  electrolytes  based  on  reverse  Monte  Carlo  modelling  of  diffraction  data. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion-conducting  glasses  are  well  known  due  to  their 
applications  in  electrochemical  devices  such  as  solid-state  micro¬ 
batteries,  fuel  cells,  chemical  sensors  and  smart  windows.  However 
mechanism  of  ionic  conductivity  in  these  materials  and  the  effect 
of  halide  dopants  is  still  not  fully  understood. 

Structural  studies  of  binary  phosphate  glasses  have  been 
reported  using  neutron  and  X-ray  diffraction  (see  reviews  by  Brow 
[1]  and  Hoppe  et  al.  [2]),  MD  simulations  [3-6]  and  spectroscopic 
techniques  [7].  These  binary  glasses  can  be  prepared  up  to  66  mol.% 
of  Li20.  Glass  forming  region  in  ternary  glass  system  LiP03-LiX 
[X  =  I,  Br,  Cl]  has  been  reported  along  with  their  ionic  conductiv¬ 
ity  (<r)  [8].  Glasses  form  up  to  30  mol.%  of  LiCl  and  33  mol.%  of  LiBr 
and  Lil.  a  increases  with  halogen  ion  size  (and  polarisability)  and 
with  increasing  LiX  content  throughout  the  glass  forming  region. 
Another  study  [9]  in  Li20-P205-LiCl  system  reports  glass  forming 
region  up  to  nearly  30  mol.%  of  LiCl  for  R  (R  =  Li20/P205)  varying 
between  0.67  to  1.38  and  up  to  35  mol.%  for  R  value  of  0.25.  The  glass 
transition  temperature  Tg  decreases  and  a  increases  as  LiCl  con¬ 
tent  increases  for  all  values  of  R  studied.  Conclusions  from  Raman 
spectroscopy  studies  of  0.3LiCl-0.7(0.58Li20-0.42P205)  glass  in 
the  same  study  indicate  that  the  addition  of  Li20  breaks  up  P-O-P 
bridges,  creating  non-bridging  oxygens  (NBOs)  which  reduce  the 
lengths  of  chains  formed  by  P-O-P  bridges,  while  the  addition  of 
LiCl  does  not  bring  about  any  such  modification  of  chain  lengths 
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in  the  glass  network,  though  it  reduces  the  Tg  of  the  parent  glass. 
In  contrast,  IR  and  XPS  studies  of  x(2LiCl)y(Li20)(100  -x-y)P205 
[x  =  0,  10,  20,  y  =  50,  60] [10]  suggest  that  the  more  open  struc¬ 
ture  of  the  halide  doped  glasses  is  achieved  by  Cl-  altering 
the  P-O-P  glass  network.  X-ray  diffraction  study  and  MD  sim¬ 
ulation  of  0.52Li20-0.48P205,  0.45Li20-0.42P205-0.13LiCl  and 
0.39Li20-0.36P205-.25LiCl  glasses  [11,12]  reports  that  the  addi¬ 
tion  of  Li20  to  P205  breaks  up  the  3D  branched  structure  creating 
more  non-bridging  oxygens  and  the  addition  of  LiCl  makes  the 
structure  less  rigid  and  creates  wider  channels  for  Li+  movement, 
thus  enhancing  a  and  decreasing  Tg.  Further  structural  reports  for 
Li20-P205-LiCl  glass  [7,13]  confirm  that  no  P-Cl  bonds  are  formed 
and  that  LiCl  dissolves  interstitially  in  LiP03  glass  network.  While 
Raman,  IR  spectroscopy  and  Reverse  Monte  Carlo  fits  of  diffrac¬ 
tion  data  accordingly  find  that  LiCl  acts  as  a  dopant  only  in  lithium 
borate  glasses  [14,15];  in  Li20-Te02-LiCl  glasses  [16,17],  the  addi¬ 
tion  of  LiCl  is  reported  to  produce  NBOs  leading  to  the  formation 
of  Te-Cl  bonds  in  the  equatorial  position  of  the  trigonal  pyramid 
Te04  unit.  7Li  MAS  NMR  studies  of  LiCl-Li20-P203  suggest  the  exis¬ 
tence  of  LiCl  aggregates  for  high  LiCl  contents  close  to  the  limit 
of  the  glass  forming  region,  and  it  is  concluded  that  these  lead 
to  the  high  mobility  of  Li+  ion  [18].  Yet  another  impedance  spec¬ 
troscopy  study  of  Li2X-B203-Si02  (X  =  0,  Cl2)  [19]  concludes  that 
the  network  dilatation  brought  about  by  the  addition  of  Cl-  has 
little  influence  on  conductivity. 

In  order  to  further  understand  the  relation  between  pathway 
volume  of  mobile  Li+  ion  and  the  influence  of  LiCl  content  in 
xLiCl-(l  -x)(0.6Li20-0.4P205)  (where  x  =  0.1,  0.15, 0.20, 0.25)  glass 
preparation  and  bond  valence  approach  was  applied  for  the  MD 
simulated  glasses  here. 
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2.  Methodology 

2 A.  Sample  preparation  and  properties  characterization 

Ternary  glasses  xLiCl-(l -x)(0.6Li20-0.4P205)  (where  x  =  0.1, 
0.15, 0.20, 0.25)  were  prepared  by  the  conventional  melt  quenching 
method  [10]  using  NH4H2P04,  Li2C03  and  LiCl.  Due  to  the  hygro¬ 
scopic  nature  of  phosphate  glasses  the  samples  were  stored  in  a 
desiccator  with  silica.  It  was  noticed  that  LiCl  partially  evaporates 
during  preparation  and  the  amount  of  Cl-  present  in  the  glass  there¬ 
fore  to  find  the  exact  composition  of  the  prepared  sample  we  used 
Scanning  Electron  Microscopy  (SEM)  with  the  aid  of  Energy  Disper¬ 
sive  X-ray  Detection  (EDX).  Table  1  shows  a  comparison  between 
the  expected  and  observed  chloride  contents  in  the  correspond¬ 
ing  glasses.  Hereafter  experimental  data  will  refer  to  the  observed 
compositions.  MD  simulations  are  conducted  for  the  expected  com- 
postions. 

The  quenched  samples  were  characterized  by  X-ray  powder 
diffractometry  using  Cu  Ka  radiation  (PANalytical  X’Pert  PRO) 
equipped  with  a  fast  linear  detector  (X’Celerator),  which  can  collect 
100  steps  at  the  same  time.  XRD  data  were  collected  in  the  20  range 
10-110°  with  a  nominal  scan  rate  of  70s  step-1  and  a  step  size  of 
0.05°  at  room  temperature.  No  sharp  Bragge  peaks  were  observed.  It 
conforms  glassy  nature  of  the  samples.  Density  of  glasses  was  mea¬ 
sured  by  Archimedes’  principle  (Density  Determination  Kit  Mettler 
ME-33360)  using  toluene  as  the  immersion  liquid.  The  glass  tran¬ 
sition  temperature,  Tg  was  determined  using  Differential  Scanning 
Calorimetry  (DSC;  TA  Instruments  2920  Modulated  DSC )  at  a  heating 
rate  of  10  K  min-1  under  constant  flow  of  high  purity  N2  gas. 

Ionic  conductivity  measurements  at  different  temperatures 
were  carried  out  by  impedance  spectroscopy  ( Schlumberger 
Solartron  SI  1260 )  in  the  frequency  range  of  1  Hz  to  1  MHz.  Commer¬ 
cially  available  silver  paint  (ion-blocking  electrode)  was  applied 
on  both  sides  of  the  glass  sample  to  achieve  a  controlled  contact 
area  for  the  ionic  conductivity  determination.  The  sample  was  kept 
in  a  Kiel-type  Electrochemical  Cell  (Ionic  Systems).  At  each  tem¬ 
perature  the  sample  was  kept  for  10  min  for  thermal  equilibration. 
The  bulk  resistance  Rb  was  determined  from  fitting  impedance  data 
to  Nyquist  plots  (- Z "  versus  Z').  The  equivalent  circuit  consists  of 
Cb,  Rb  and  CPE,  where  Cb  is  the  geometric  capacitor  between  the 
electrodes,  and  CPE  is  the  constant  phase  element  due  to  the  polar¬ 
isation  distribution  at  the  interface  between  blocking  electrodes 
and  glass. 


r  =  1.47  A 


Table  2 

Optimized  two-body  potential  parameters  for  LiCl  doped  phosphate  glasses 


1500, 1000  and  700  K,  equilibration  time  was  300  ps  followed  by  a 
NVT  simulation  (Hoover  thermostat)  for  200  ps.  At  500, 400, 300  K 
the  equilibration  time  was  increased  to  600  ps  and  the  production 
run  after  equilibration  was  also  increased  to  400  ps.  Thus  the  calcu¬ 
lation  of  mean  square  displacement  (MSD)  for  all  atoms  was  based 
on  the  final  200  ps  at  1500, 1000, 700  K,  and  the  final  400  ps  for  the 
temperatures  500,  400  and  300  K  below  the  glass  transition  tem¬ 
perature.  Besides  the  Coulomb  interactions  the  forcefield  involves 
a  Born  repulsion  term  (Buckingham  potential  with  Cy  =  0)  and  a 
Vessal-type  three-body  term  as  follows: 


ZZe^  /  -n, 

4>ij  =  — ~ —  +  Ay  exp  (  — 
Pij  0 


2.2.  Computer  simulations 

Constant  volume  (NVT)  MD  simulations  have  been  performed 
on  structure  models  of  (1 -x)(0.6Li20-0.4P205)-xLiCl  glasses 
(where  x  =  0.10,  0.15,  0.20,  0.25).  The  volume  of  the  simulation  box 
is  derived  from  our  density  determinations  (Fig.  2).  The  optimized 
potential  parameters  shown  in  Tables  2  and  3  have  been  derived  in 
the  frame  of  this  project  by  fitting  to  experimentally  known  bond 
lengths  and  coordination  numbers. 

Starting  at  a  temperature  of  1500  K,  the  system  containing 
1975-2170  atoms  was  cooled  by  scaling  velocities  to  1000,  700, 
500,  400  K  and  finally  cooled  to  300  K  with  a  time  step  of  2  fs.  At 

Table  1 

Chloride  content  (wt%)  in  the  studied  (1  -x)(0.6Li20-0.4P205)-xLiCl  glasses 


+  8C0^{[(0o-7r)2-(^-7r)2]2}eXP 

Here  (j)y  is  the  potential  function  between  particles  i  and  j,  Zz  is 
the  effective  atomic  charge  of  atom  i,  00  is  the  equilibrium  angle 
subtended  by  ry  and  rik.  Ay,  py0,  ky  and  py  are  forcefield  param¬ 
eters.  The  ionic  charges  used  in  MD  simulations  for  Li,  Si,  P,  O,  Cl 
are  +0.6e,  +2.4e,  +3e,  -1.2e,  -0.6e,  respectively.  Fig.  1  shows  the 
pair  correlation  function  (PCF)  and  Running  coordination  (RCN) 
for  x  =  0.25  glass,  indicating  the  presence  of  P04  tetrahedra.  Bond 

Table  3 

Potential  parameters  for  the  three-body  Vessal  term  in  the  forcefield  for  LiCl  doped 
phosphate  glasses 


Interaction  ( i-j-k )  O-P-O 

k  (eV)  0.0129 

0Q  103° 

Pij  (A)  0.2252 

pik  (A)  0.1191 

rc*  (A)  1.90 


rc*:  cut-off  in  r,j  and  rjk. 


rJ_  _rJk_ 
Pij  Pjk 
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valence  analysis  is  applied  to  the  final  configuration  obtained  at 
300  K. 

2.3.  Bond  valence  approach 

Empirical  relationships  between  bond  length  and  bond  valence 
(BV)  are  widely  used  to  identify  plausible  equilibrium  sites  for 
an  atom  in  a  crystal  structure  as  sites  where  the  BV  sum  of  the 
atom  matches  its  oxidation  sate  [20].  A  modification  of  the  BV 
approach  that  improves  the  assessment  of  non-equilibrium  site 
energies  by  a  systematic  adjustment  of  BV  parameters  to  the  bond 
softness  [21-26]  extends  the  application  range  of  this  simple  tool 
to  structure-property  relationship  studies  in  solid  electrolytes.  For 
that  purpose,  ‘accessible’  sites  for  a  mobile  ion,  A,  in  a  local  struc¬ 
ture  model  are  identified  using  empirical  relationships  between  the 
bond  length,  R  and  a  so-called  bond  valence,  s^_x, 


sA-x  =  exp 


r(Ko-K)l 

b 


Fig.  2.  Composition  dependence  of  the  mass  density  (solid  symbols)  and  number 
density  (open  symbols)  in  the  glasses  of  (1  -x)(0.6Li20-0.4P205)-xLiCl.  The  lines 
are  drawn  as  a  guide  to  the  eye. 


as  sites  where  the  mismatch  of  the  bond  valence  sum,  V(A) 

I A  17(A) |  =  |^>_x  -  Vid(A)|  +  y>_x 
x  x 

over  the  sA_x  from  all  adjacent  counterions,  X  approaches  the  ideal 
valence,  Vid(A)  (which  equals  its  oxidation  state).  To  enhance  the 
chemical  plausibility  of  BV  mismatch  ‘energy  landscapes’,  |AV(A)| 
contains  a  penalty  function,  pA_x,  that  discriminates  against  sites 
where  a  matching  V(A)  is  achieved  by  unfavourable  strongly  asym¬ 
metric  coordinations  [24]  and  is  complemented  by  hard  minimum 
distance  and  soft  coordination  number  constraints.  In  contrast  to 
most  conventional  BV  parameter  sets  our  softBV  parameters  are 
based  on  the  assumption  that  not  only  the  counterions  of  the  first 
coordination  shell  but  all  counterions  up  to  a  cut-off  radius  of  4-8  A 
(depending  on  sizes  and  polarisabilities  of  the  ions  involved)  con¬ 
tribute  to  V(A).  This  is  indispensable  for  modelling  ion  transport 
pathways  as  it  avoids  artefacts  in  the  BV  sum  variation,  when  an  ion 
moves  across  the  border  of  its  coordination  shell.  BV  ion  transport 
pathway  models  are  based  on  the  supposition  that  paths  between 
accessible  sites,  along  which  |  AV(A)|  remains  sufficiently  low  rep¬ 
resent  probable  ion  transport  paths.  Isosurfaces  of  fixed  maximum 
BV  sum  mismatch  |  AV(A)|  for  a  certain  ion  type  A  will  thus  enclose 
regions  that  an  A  ion  in  the  pathway  may  reach  with  a  certain 
activation  energy,  Fa.  Low  BV  mismatch  regions  that  include  both 
occupied  and  vacant  sites,  therefore,  enable  local  jumps  of  A.  A 
long-range  transport  requires  pathways  that  extend  through  the 
whole  structure.  For  periodic  models  this  is  equivalent  to  span¬ 
ning  the  model  box.  Determining  the  ‘pathway  volume’  for  a  given 
maximum  value  of  |AV(A)|  yields  a  simple  and  reliable  way  of 
characterizing  ion  transport  pathways  in  local  structure  models  of 
disordered  systems  [23,25].  To  evaluate  this  pathway  volume,  V(A), 
has  to  be  calculated  for  a  hypothetical  A+  ion  at  each  point  in  the 
structure.  In  practice,  the  structure  model  is  divided  into  a  prim¬ 
itive  grid  containing  ~3.3  million  cubic  volume  elements  with  a 
size  of  ca.  (0.2  A)3,  i.e.  about  1600  volume  elements  per  atom.  A 
volume  element  is  considered  to  belong  to  the  conduction  path¬ 
way  if  |V(A)  —  V(A)ideai|  is  below  a  fixed  threshold  value  or  if  the 
sign  of  (V- Videal)  alters  within  the  volume  element  (the  second 
condition  cushions  effects  of  the  limited  grid  resolution).  Accessi¬ 
ble  volume  elements  that  share  common  faces  or  edges  belong  to 
the  same  ‘pathway  cluster’.  It  is  assumed  that  these  pathway  clus¬ 
ters  contribute  to  dc  conductivity  only  if  they  percolate  through  the 
structure  model.  The  volume  fraction,  F,  of  the  percolating  path¬ 
ways  turns  out  to  be  a  suitable  measure  of  Fa  and  thereby  (in  the 
case  of  glasses)  also  of  the  room  temperature  conductivity.  The 


remaining  restricted  pathway  clusters  are  considered  to  contribute 
only  to  ac  conductivity.  As  mentioned  above,  the  threshold  value 
for  |  AV(A)|  has  to  be  chosen  somewhat  arbitrarily.  When  compar¬ 
ing  pathways  for  different  types  of  mobile  ions,  we  find  that  the 
most  suitable  value  will,  as  discussed  further  below,  depend  on  the 
mass,  mA ,  of  the  mobile  ion. 

3.  Results  and  discussions 

Both  measured  density  and  number  density  of  glasses  are 
decreasing  with  increase  of  LiCl  content  (Fig.  2).  This  indicates 
that  the  glass  network  is  becoming  more  open  which  leads  to  the 
decrease  of  Tg  (Fig.  3)  and  an  increase  of  Li+  ion  mobility.  Complex 
impedance  measurements  of  the  glasses  were  carried  out  from  27 
to  100  °C.  The  Nyquist  plots  were  fitted  by  using  an  equivalent  cir¬ 
cuit  model  of  one  RC  couple  (representing  the  bulk  resistance  and 
geometric  capacitance)  in  series  with  one  CPE  (representing  polar¬ 
isation  at  the  electrode:  electrolyte  interface)  (Fig.  4).  A  series  of 
Nyquist  plots  of  0.54Li20-0.36P205-0.10LiCl  glass  at  different  tem¬ 
peratures  is  shown  in  Figs.  4  and  5  shows  the  Arrhenius  behaviour  of 
dc  ionic  conductivity  as  expected.  Table  4  shows  the  experimental 
and  MD  ionic  conductivity  and  at  27  °C.  Comparison  of  ionic  con¬ 
ductivity  of  both  experimental  and  simulations  is  shown  in  Fig.  6. 
MD  values  are  of  the  same  order  of  the  experimental  values.  The 


Fig.  3.  Composition  dependence  of  the  glass  transition  temperature  (Tg)  in  the 
glasses  of  (1  -x)(0.6Li20-0.4P205)-xLiCl. 
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Fig.  4.  Nyquist  plots  of  0.54Li20-0.36P205-0.10LiCl  glass  at  different  temperatures 
and  equivalent  circuit  used  for  impedance  fitting. 


Fig.  5.  Arrhenius  plots  of  the  dc  conductivities  obtained  from  impedance  spec¬ 
troscopy  for  different  glasses. 

trend  of  increment  of  conductivity  is  same  in  both  cases.  The  slope 
of  Arrhenius  plots  decreases  with  increase  of  LiCl,  consequently 
activation  energy  (Ea)  decreases  form  0.494  eV  (47.6  kj  mol-1 )  to 
x  =  0.25  (Fig.  7). 

Fig.  8  shows  the  presence  of  P04  tetrahedra  in  MD  generated 
configuration  of  0.45Li20-0.30P205-0.25LiCl  glass  at  300 K.  Final 
equilibrated  MD  configurations  glasses  were  used  for  bond  valence 
analysis.  The  threshold  value  for  |  A  V(A)|  is  chosen  to  be  0.10  valence 
units  to  allow  for  a  comparison  with  existing  reference  data  for  a 
wide  range  of  other  oxide  glasses.  Using  the  bond  valence  approach 
it  can  be  investigated  whether  the  Li+  ions  form  LiCl  clusters  in 
the  investigated  glass  structures,  as  had  been  discussed  previ¬ 
ously  in  the  literature  for  related  systems  [24].  It  is  found  that  this 
hypothesis  can  be  ruled  out  for  the  system  under  study:  the  dis¬ 
tribution  of  relative  contributions  of  Li-Cl  bonds  to  the  Lithium 
BV  sums  in  these  glasses  (cf.  Fig.  9)  demonstrates  that,  as  to  be 


Fig.  6.  Comparison  of  ionic  conductivities  for  ( 1  -  x)(0.6Li20-0.4P205  )-xLiCl  glasses 
with  composition  at  27  °C.  The  lines  are  linear  fits  of  data. 


Fig.  7.  Composition  dependence  of  activation  energy  (£a)  in  the  glass  series 
(1  -x)(0.6Li20-0.4P205)-xLiCl. 

expected  for  a  modified  random  network  glass,  only  a  small  frac¬ 
tion  of  the  Li  atoms  are  bonded  predominantly  by  chloride  ions 
(<2%  for  all  glasses).  For  none  of  the  Li  atoms  in  these  glasses  the 
coordination  polyhedron  exists  exclusively  of  Cl-  ions.  More  than 
95%  oxide  contribution  occurs  for  20%  (x  =  0.25;  0:C1  ratio  7.8:1)  to 
64%  (x  =  0.10,  0:C1  ratio  23.4:1)  of  the  Li  atoms.  This  also  clarifies 
that  any  lithium  ion  transport  pathway  has  to  run  along  sites  with 
mixed  oxide/halide  or  exclusive  oxide  coordination.  Fig.  10  shows 
the  local  environment  of  the  Li+  ion  sites  located  in  the  conduc¬ 
tion  pathways  of  glasses  under  present  investigation.  The  density  of 
pathways  slightly  increases  with  increasing  chloride  content  from 
x  =  0.10  to  0.25  (the  volume  fraction  of  percolating  pathways  (F) 
increases  from  5.4%  to  5.8%),  in  line  with  the  increase  of  ionic  con¬ 
ductivity  in  these  glasses,  except  at  x  =  0.15.  Figs.  11  and  12  relate 
the  volume  fraction  of  the  percolating  ion  conduction  pathways  in 
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Fig.  8.  Structure  of  0.45Li20-0.30P205-0.25LiCl  glass  at  300  K.  Oxide  atoms  (red 
spheres)  are  around  P  (shown  as  tetrahedra),  Li:  blue  spheres,  Cl:  green  spheres. 
Only  1/4  of  the  structure  model  is  shown  along  z  (perpendicular  to  the  plane  of 
view)  to  reduce  overlap.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


Fig.  9.  Distribution  of  relative  contributions  of  Li-Cl  bonds  to  the  BV  sums  of  Li+  ions 
in  the  glasses  (1  -x)(0.6Li20-0.4P205)-xLiCl  for  x  =  0.10-0.25.  For  the  case  x  =  0.25 
only  10%  of  the  Li  experience  a  BV  contribution  exceeding  34%. 


Fig.  11.  Variation  of  pathway  volume  fraction  of  Li+  ion  with  experimental 
room  temperature  ionic  conductivity.  Data  in  solid  symbol  refer  to  data  from 
RMC  models  [23].  Data  in  open  circles  refer  to  MD  simulated  structures  of 
(1  -x)(0.6Li20-0.4P205)-xLiCl  glasses  (where  x  =  0.10,  0.15,  0.20  and  0.25). 


Fig.  12.  Variation  of  pathway  volume  fraction  of  Li+  ion  with  activation  energy.  Data 
in  solid  symbol  refer  to  data  from  RMC  models  [23].  Data  in  open  circles  refer  to  MD 
simulated  structures  of(l  -xX0.6Li2O-0.4P2Os)-xLiCl  glasses  (where  x  =  0.10,  0.15, 
0.20  and  0.25). 


Fig.  10.  Slices  of  the  isosurfaces  of  constant  lithium  bond  valence  sum  mismatch  |  AV(Li)|  in  the  glasses  (1  -x)(0.6Li20-0.4P205)-xLiCl  for  x  =  0.15,  0.20,  0.25  (left  to  right)  at 
300  K,  projected  along  the  z-axis. 
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Fig.  13.  The  local  Li+  ion  transport  pathway  dimension,  Dm(r)  vs.  radius  (r)  for 
(1  -x)(0.6Li20-0.4P205)-xLiCl  glasses.  Inside  graph  indicates  the  variation  of  local 
dimensional  minima  with  respect  to  LiCl  variation. 

the  investigated  glasses  to  their  experimentally  observed  absolute 
values  of  the  ionic  conductivity  and  activation  energy  (Ea).  Results 
of  the  present  MD  study  are  compared  with  the  previously  reported 
RMC  values.  As  seen  from  Figs.  11  and  12,  the  increase  in  the  vol¬ 
ume  fraction  when  we  move  from  the  glass  x  =  0.1  to  0.25  leads  to 
an  increase  in  the  ionic  conductivity  and  a  corresponding  decrease 
in  the  activation  energy,  except  for  x  =  0.15.  Using  the  concept  of 
mass-radius  dimension  [14],  we  can  further  investigate  details  of 
the  pathway  structure.  As  seen  from  Fig.  13,  the  minimum  value 
of  local  pathway  dimension  increases  with  the  LiCl  content,  again 
except  for  x  =  0.15.  This  minimum  value  of  local  dimension  may  be 
largely  thought  of  as  characterizing  the  width  of  the  Li  ion  transport 
pathway  at  the  bottleneck  of  elementary  transport  steps.  Thus  the 
increase  of  the  local  pathway  dimension  as  a  result  of  doping  by 
the  more  polarisable  Cl-  ion  indicates  that  channels  accessible  for 
Li+  ion  movement  are  widened  by  the  doping.  This  contributes  to 
the  decrease  of  activation  energy  and  increase  of  ionic  conductivity 
and  with  rising  LiCl  content. 

4.  Concluding  remarks 

Optimized  potential  parameters  reproduce  the  reported  bond 
lengths  and  coordinations.  Moreover,  experimental  and  simulated 


values  of  the  ionic  conductivity  are  well  in  agreement.  The  effect  of 
LiCl  addition  on  the  pathways  for  the  mobility  of  Li+  ions  is  quan¬ 
tified  for  the  phosphate  glasses  xLiCl-(l  -x)  (0.6Li20-0.4P205). 
Pathway  volume  fraction  for  Li  ion  transport  and  experimental 
conductivities  accordingly  exhibit  an  unusually  small  effect  of 
halide  doping  on  the  conductivity.  For  the  investigated  composi¬ 
tions  ion  transport  cannot  be  related  to  (previously  postulated) 
LiCl  aggregates  as  such  aggregates  do  not  exist— at  least  in  the 
simulated  structures.  Nearly  all  Li  have  a  mixed  oxide  chloride 
coordination.  The  minimum  value  of  local  pathway  dimension 
increases  with  the  LiCl  content,  except  for  x  =  0.15.  The  higher 
minimum  value  of  local  dimension  means  the  Li  ion  transport  path¬ 
way  is  more  open,  i.e.  channels  for  the  Li  ion  movement  become 
wider. 
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